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CN ' Abstract 

D | This paper presents the analysis of Doppler p- modes of the F9 V star /3 Virginis ob- 

tained with the spectrograph Coralie in March 2003. Eleven nights of observations 
have made it possible to collect 1293 radial velocity measurements with a standard 
deviation of about 2.2 m s _1 . The power spectrum of the high precision velocity time 
£> ' series clearly presents several identifiable peaks between 0.7 and 2.4 mHz showing 

regularity with a large and small spacings of Av = 72.1 fiHz and 5uq2 = 6.3 /iHz re- 
spectively. Thirty-one individual modes have been identified with amplitudes in the 

CN ' range 23 to 46 cms -1 , i.e. with a signal to noise between 3 and 6. 
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1 Introduction 



The lack of observational constraints leads to serious uncertainties in the mod- 
eling of stellar interiors. The measurement and characterization of oscillation 
modes is an ideal tool to test models of stellar inner structure and theories 
of stellar evolution. Indeed, the five-minute oscillations in the Sun have led 
to a wealth of information about the solar interior. These results stimulated 
various attempts to detect a similar signal on other solar-like stars. Solar-like 
oscillation modes generate periodic motions of the stellar surface with periods 
in the range of 3-60 minutes but with extremely small amplitudes. Essen- 
tially, two methods exist to detect such a motion: photometry and Doppler 
spectroscopy. In photometry, the oscillation amplitudes of solar-like stars are 
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within 2-30 ppm, while they are in the range of 10-150 cms -1 in radial ve- 
locity measurements. Photometric measurements made from the ground are 
strongly limited by scintillation noise. To reach the needed accuracy requires 
observations made from space. In contrast, Doppler ground-based measure- 
ments have recently shown their ability to detect oscillation modes in solar-like 
stars. These past years, the stabilized spectrographs developed for extra-solar 
planet detection achieved accuracies needed for solar-like oscillation detection 
by means of radial velocity measurements (Carrier et al. 2003). 

A primary target for the search for p-mode oscillations is the bright dwarf 
F9 Virginis (HR 4540, HD 102870, m v = 3.61). Scaling from the solar case, 
p-modes are expected near z/ max = 1.4mHz with a large frequency spacing of 
about Az/ = 71 /iHz, and a maximal amplitude of A osc = 65cms _1 (Kjeldsen 
& Bedding 1995). 

In this paper, we report Doppler observations of (3 Vir made with the Coralie 
spectrograph, well known for the characterization of p-modes on the a Cen 
system (Bouchy & Carrier 2002, Carrier & Bourban 2003). These measure- 
ments enable the identification of thirty-one individual mode frequencies. 



2 Observations and data reduction 

(3 Vir was observed over a campaign of eleven nights (2003 February 28 - 
March 10) with Coralie, the high-resolution fiber-fed echelle spectrograph 
mounted on the 1.2-m Swiss telescope at La Silla (ESO, Chile). During the 
stellar exposures, the spectrum of a thorium lamp carried by a second fiber is 
simultaneously recorded in order to monitor the spectrograph's stability and 
thus to obtain high-precision velocity measurements. The wavelength coverage 
of the spectra is 3875-6820 A, recorded on 68 orders. The dead-time between 
two exposures was improved from the usual 125 s to 85 s, by archiving the 
image during the following exposure. Exposure times of 120 s, thus cycles of 
205 s, allowed us to obtain 1293 spectra, with a typical signal-to-noise ratio 
(S/N) in the range of 100-140 at 550 nm. 

Radial velocities are computed for each night relative to the highest S/N spec- 
trum obtained in the middle of the night by the use of the optimum-weight 
procedure (Connes 1985, Carrier et al. 2001). This method requires a Doppler 
shift that remains small compared to the line-width (smaller than 100ms -1 ). 
Since the Earth's motion can introduce a Doppler shift larger than 700ms _1 
during a whole night, each spectrum is first corrected for the Earth's motion 
before deriving the radial velocities. Subsequently, the mean for each night is 
subtracted. The rms scatter of the time series is 2.24 ms" 1 (see Fig. 1 and Ta- 
ble 1), which can be directly compared with the fundamental uncertainty due 
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Fig. 1. Radial- velocity measurements of (3 Vir. The dispersion (which includes the 
noise and the oscillations) reaches 2.24 ms" 1 . The upper pannel represents the pho- 
ton noise uncertainties. 



Table 1 

Distribution and dispersion of Doppler measurements 



Date 


Nb spectra 


Nb hours 


cr (ms r ) 


2003/02/28 


69 


4.67 


1.91 


2003/03/01 


91 


6.06 


2.35 


2003/03/02 


118 


8.21 


2.01 


2003/03/03 


137 


8.55 


2.25 


2003/03/04 


135 


7.95 


2.22 


2003/03/05 


142 


8.37 


2.33 


2003/03/06 


111 


6.75 


2.29 


2003/03/07 


132 


8.24 


2.60 


2003/03/08 


76 


9.13 


2.54 


2003/03/09 


146 


8.64 


2.12 


2003/03/10 


136 


7.92 


2.03 
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Fig. 2. Power spectrum of the Coralie radial velocity measurements of j3 Vir. The 
window function is shown in the inset. 

to photon noise. The uncertainties coming from the thorium spectrum used 
in the instrumental tracking are quite stable with a value of about 0.34 ms 4 . 
The quadratic sum of the stellar spectrum photon noise and this instrumental 
photon noise varies between 1.2 and 2.3 ms -1 . 



3 Power spectrum analysis 



3.1 Coralie measurements 



In order to compute the power spectrum of the velocity time series, we use the 
Lomb-Scargle modified algorithm (Lomb 1976, Scargle 1982) with a weight 
being assigned to each point according its uncertainty estimate. The time 
scale gives a formal resolution of 1.14 //Hz. The Nyquist frequency (cut-off) is 
determined by the measurement cycle of 205 s and has a value of ~ 2.4 mHz. 
The resulting periodogram, shown in Fig. 2, exhibits a series of peaks centred 
at 1.4 mHz, exactly where the solar-like oscillations for this star are expected. 
Typically for such a power spectrum, the noise has two components: 

• Towards the lowest frequencies, the power should scale inversely with fre- 
quency squared, as expected for instrumental instabilities. However, the 
computation of the radial velocities introduces a high pass filter. Indeed, 
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Fig. 3. Power spectrum of the Harps radial velocity measurements of (5 Vir (see 
Sect 3.2). The daily aliases are too high to allow any mode determination, however, 
we can see that p-modes are present in the range 0.7-2.4mHz. 

the radial velocities were computed relative to one reference for each night 
and the average radial velocities of the night fixed to zero (see Sect. 2). This 
results in an attenuation of the very low frequencies which can be seen on 



• At high frequencies it is flat, indicative of the Poisson statistics of photon 
noise. As p-modes are present until a frequency of 2.4 mHz (see Sect. 3.2), 
we have to compute this noise at lower frequency than the oscillations. The 
noise was determinated per interval: it is decreasing until 530 /iHz (instru- 
mental instabilities) and is increasing again from 690 /xHz due to the oscil- 
lation modes. We thus calculated the flat noise in the interval 530 - 690 /iHz 
which seems constant with a value of 0.0074 m 2 s~ 2 , namely 7.6 cm s -1 in 
amplitude. With 1293 measurements, this high frequency noise corresponds 



The power spectrum presents an excess in the range 0.7-2.4 mHz (see also 
Sect. 3.2). Note that the filtering induced by the radial velocities computation 
does not influence the frequency of the peaks in the range 0.7-2.4 mHz, but 
could slightly change their amplitudes. The amplitude of the strongest peaks 
reaches 46 cms -1 , corresponding to a signal to noise of 6 (in the amplitude 
spectrum) . 

3.2 Harps measurements 

Some measurements were obtained with the spectrograph Harps installed on 
the 3.6-m telescope at La Silla Observatory (ESO, Chile) (Pepe et al. 2002) in 
July 2004 in order to better determine the range of p-mode frequencies. 171 
measurements were collected during 5 nights and the resulting power spectrum 
is shown in Fig. 3. Due to observations of only 40 minutes per night, these 
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Fig. 4. Autocorrelation of the power spectrum with a threshold of 0.07 m 2 s _1 . The 
gray line corresponds to the same autocorrelation but all peaks greater than the 
threshold have the same amplitude (not to favour only the highest peaks). The 
large spacing is estimated to be about 72 fiUz. 

measurements cannot allow us to determine mode frequencies: the five first 
daily aliases (until 5 x 11.57/iHz on each side) have an amplitude greater than 
90% of the main mode: it is impossible to safely choose the right frequencies. 
However, due to a quicker observation cycle and thus to a higher Nyquist 
frequency, we can remark that oscillation modes are present in the range 0.7- 
2.4 mHz. 



3.3 Search for a comb-like pattern 



In solar-like stars, p-mode oscillations of low-degree are expected to produce a 
characteristic comb-like structure in the power spectrum with mode frequen- 
cies v n> i reasonably well approximated by the asymptotic relation (Tassoul 
1980): ' 



is n>l nAv(n + - + e)-l(l + l)D . (1) 

Here, D , which is equal to ^5v 02 if the asymptotic relation holds exactly, 
and e are sensitive to the sound speed near the core and to the surface layers 
respectively. The quantum numbers n and I correspond to the radial order 
and the angular degree of the modes, and Av and to the large and small 
spacings. To fit to this relation, an autocorrelation of the power spectrum 
is calculated and presented in Fig. 4. The threshold of the power spectrum 
(all frequencies with a smaller amplitude than this value has its amplitude 
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fixed to zero) is chosen to 0.07m 2 s _1 corresponding to a signal to noise of 
3.5 in amplitude, this process can diminish the noise. Each peak of the au- 
tocorrelation corresponds to a structure present in the power spectrum. The 
two strong peaks at low frequency near 11.6 and 23.1 /iHz correspond to the 
daily aliases. The strongest other peaks are situated near 60 ± 11.57 /iHz and 
130 ± 11.57 /iHz which correspond respectively to one and two times the large 
spacing. We can thus deduce that the large spacing is near 72 /^Hz, as peaks 
at 2 times a separation of 60 or 80 /iHz are not present, this is in accordance 
with theoretical predictions (see Sect. 1). Three peaks coexist near 72 /iHz, 
namely 70.5, 72 and 74 /iHz, but the most probable large spacing is situated 
at 72 /iHz: a high peak has the value of 2 times this spacing whereas only a 
low peak is present near 141 /iHz, and no peak at all near 148 /^Hz. This large 
spacing will be in any case refined during the following step. 



3.4 Mode identification 

The frequencies were extracted using an iterative algorithm, which identifies 
the highest peak between 700 and 2400 /iHz and subtracts it from the time 
series. Note that because of the stochastic nature of solar-like oscillations, a 
timestring of radial velocities cannot be expected to be perfectly reproduced 
by a sum of sinusoidal terms. Therefore, using an iterative clean algorithm to 
extract the frequencies can add additional peaks with small amplitudes due 
to the finite lifetimes of the modes that we do not know. Nevertheless, the 
iterative algorithm ensures that one peak and its aliases with an amplitude 
above a given threshold is only extracted once. To avoid extracting artificial 
peaks with small amplitudes added by the iterative algorithm, the choice of 
this threshold is important. In the case of (5 Vir, we iterated the process until 
all peaks with an amplitude higher than 3 a in the amplitude spectrum were 
removed (see Fig. 5). Peaks with amplitudes below the 3a threshold were not 
considered since they were too strongly influenced by noise and by interactions 
between noise and daily aliases. This threshold, which ensures that the selected 
peaks have only a small chance to be due to noise, gave a total of twenty-three 
frequencies (see Table 2). Because of the daily alias of 11. 57 /iHz introduced 
by the monosite observations, we cannot know a priori whether the frequency 
selected by the algorithm is the right one or an alias. We thus considered that 
the frequencies could be shifted by ±11.57 /iHz, and made echelle diagrams for 
different large spacings near 72 /iHz until every frequency could be identified 
as an £ = 0, £ = 1 or £ = 2 mode. In this way, we found an averaged large 
spacing of 72.1 /iHz. 

To investigate how many peaks are expected to be due to noise, we have con- 
ducted simulations in which we analyzed noise spectra containing no signal. 
For this purpose, a velocity time series is built, using the observational time 
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Table 2 

Identification of extracted frequencies. The number of peaks due to noise is in agree- 
ment with the simulations described in Sect. 3.4, which predict 4.9 ± 4.3 noise peaks. 



Frequency 
[/xHz] 


Mode ID 


S/N 


von 1 i 1 i e 
t zU.l + 11. b = 


7Q1 1 

- 1 ol. / 


£ = 1 


3.2 


Q1 A R 1 1 

814. b — 11. b = 


qhq n 
= 8U0.U 


£ = 1 


3.0 


non n 1 1 e 

yzu.y — ii. b = 


= yuy.o 


^ = 


3.3 




non e 

982.0 


= 


4.2 


mo o i 1 1 e 
lllo.o + 11.0 = 


1 1 i/i n 

iiz4.y 


= 


4.1 


lice q ii e 

llbo.o — 11. — 


1 1 cq 7 


£ = 1 


3.8 


1010 l 11 c 

1Z1Z.0 + 11. — 


1 OO/I o 
1ZZ4.Z 


£ = 1 


3.1 




1 OKC 7 


£ = 2 


3.3 




1 Ofifi 1 

Izbb. 1 


£ = 


4.1 




1 om T 

Izyl. / 


noise 


4.1 




1 one 7 

izyo. / 


£ = i 


4.6 


i Qon 1 lie 

lozy.i + 11.0 = 


1340. / 


= 


3.9 


1 QOO O lie 
lOOZ.Z — 11.0 = 


1 Q^n e 
13 lU.b 


£ = 1 


4.6 


1 QH1 Q I 1 1 e 

loyi.o + ii.o = 


1 /i nq o 
14Uo.z 


£ = 2 


3.0 


1/1010 1 1 e 
14Z1.Z — 11.0 = 


1 /l nn e 
14(jy.b 


£ = 


6.1 


1 Jon n i lie 

143U. y + 11. D = 


144Z.0 


£ = 1 


3.9 




144/ .y 


noise 


3.3 




14 1<S.0 


£ = 2 


4.2 




1 /i vn n 

14/y.y 


£ = 2 


4.5 


1 C^07 Q 1 1 e 
lDZV.8 — 11.0 = 


lolb.z 


£ = 1 


4.5 




i e 
looo.b 


1 = 2 


5.2 


i enc v i 1 e 
ioy8. / — 11.0 = 


1 rov 1 

108/ .1 


1=1 


3.3 




1 eo/1 1 
lbz4. 1 


1 = 2 


4.0 




1697.8 


1 = 2 


3.0 


1691.1 + 11.6 = 


1702.7 


I = 


3.1 


1754.9 + 11.6 = 


1766.5 


^ = 2 


3.1 


1818.9 - 11.6 = 


1807.3 


1 


3.9 


1835.0 + 11.6 = 


1846.6 


£ = 


3.5 




1875.2 


£= 1 


3.6 


2147.0 - 11.6 = 


2135.4 


^ = 


4.1 




2168.4 


1 


3.2 
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0.15 
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1000 



1500 
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2000 



Fig. 5. Top: Original power spectrum of (3 Vir. Bottom: Cleaned power spectrum: 
all peaks listed in Table 2 have been removed. The dot-dashed, dotted and dashed 
lines indicate an amplitude of 5 a, 4<r and 3 c, respectively. Numerous peaks are 
still present below 3 a, since no peaks have been cleaned below this threshold. These 
peaks can be due to p-mode oscillations and noise or have artificially been added 
by the extraction algorithm due to the finite lifetimes of the modes 

sampling and radial velocities randomly drawn by assuming a Gaussian noise 
(Monte-Carlo simulations). The amplitude spectrum of this series is then cal- 
culated and peaks with amplitude greater than 3, 4 and 5 a are counted; note 
that a peak and its aliases are only counted once. The whole procedure is 
repeated 1000 times to ensure the stability of the results. In this way, we find 
that the number of peaks due to noise with an amplitude larger than 3 a is 
4.9 ±4.3 in the range 0.7-2.4 mHz, for 4a, the number of peaks due to noise 
varies between and 3 with a mean value of 0.0 and a standard deviation of 
0.5. No peaks due to noise are expected with an amplitude larger than 5 a. 



The echelle diagram showing the thirty-one identified modes is shown in Fig. 7. 
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Table 3 

Oscillation frequencies (in ^Hz) of (5 Vir. The frequency resolution of the time series 
isl.l4/xHz. 







I = 


£ = 1 


I = 2 


n = 


8 




731.7 




n = 


9 




803.0 




n = 


10 








n = 


11 


909.3 






n = 


12 


982.6 






n = 


13 








n = 


14 


1124.9 


1153.7 




n = 


15 




1224.2 


1258.7 


n = 


16 


1266.1 


1298.7 




n = 


17 


1340.7 


1370.6 


1403.2 


n = 


18 


1409.6 


1442.5 


1478.5 / 1479.9 


n = 


19 




1516.2 


1553.6 


n = 


20 




1587.1 


1624.1 


n = 


21 






1697.8 


n = 


22 


1702.7 




1766.5 


n = 


23 




1807.3 




n = 


24 


1846.6 


1875.2 




n = 


25 








n = 


26 








n = 


27 








n = 


28 


2135.4 


2168.4 




n = 


29 




2238.2 




n = 


30 








n = 


31 


2352.5 






Ave 




72.18 (0.09) 


71.91 (0.15) 


72.82 (0.30) 
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Fig. 6. Power spectrum of f3 Vir with the thirty-one extracted frequencies indicated 
by dashed lines. The identification of each extracted frequency is given in Table 3. 

The frequencies of the modes are shown in Fig. 6 and are given in Table 3, with 
radial order of each oscillation mode deduced from the asymptotic relation (see 
Eq. 1) assuming that the parameter e is near the solar value (e ~ 1.5). The 
average large spacing and the parameters D Q and e are thus deduced from a 
least-squares fit of this equation with the frequencies of Table 3. The average 
small spacing is here defined as D Q x 6 (the perfect asymptotic case, with 
constant large and small separations): 

Az/ = 72.07 ± 0.10 nHz, e = 1.59 ± 0.03 , 

D = 1.05 ± 0.24 nHz, (Su 02 ) = 6.28 ± 1.42 fiHz. 



The small spacing appears to decrease from 7.4 to 4.9 /iHz between 1260 and 
1700 yuHz. (3 Vir has a projected rotational velocity of about v sini ~4.3kms _1 
(Glebocki & Stawikowski 2000). We can estimate its radius to R = 1.66R Q 
from L = 3.51L Q (Hipparcos data) and T e fj = 6140°K (Gray et al. 2001). As- 
suming an uniform rotation, we determine the corresponding splitting of the 
modes to be at least 0.56 /iHz (for sini = l). This splitting can increase the 
uncertainty of the mode frequencies. £ = 1 modes can thus be shifted by at 
least ±0.56 /iHz, this effect is worse on £ = 2 modes (twice this value). Indeed, 
the dispersions of £ = 2 and £ = 1 modes are greater than for £ = modes 
(see Fig. 7 and Table 3 with the larger errors on Az/ 2 and Aui). Moreover, 
two high £ = 2 modes were determined with a same radial order (1478.5 and 
1479.9 /iHz), this could be explained by such a splitting, implying an angle % 
near 53°, however, with a high uncertainty. 

Both £ = 1 modes at 731.7 and 803.0 /iHz could be shifted by 11.57/xHz. They 
were fixed to their value in order to follow the slight curvature present for 
£ = modes at low frequency. They have to be taken with some caution. 
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Fig. 7. Echelle diagram of identified modes with a large separation of 72.1 /iHz. The 
modes 1=2 (■), £=0 (•), and 1=1 (A) are represented with a size proportional to 
their amplitude. 
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3.5 Oscillation amplitudes 



Concerning the amplitudes of the modes, theoretical computations predict 
oscillation amplitudes near 50 cms" 1 for a 1.25 M & star like (3 Vir, with mode 
lifetimes of the order of several days to several tenth of days. (Houdek et 
al. 1999). The amplitudes of the highest modes, in the range 32-46 cms" 1 , are 
then slightly lower than expected. The observations indicate that oscillation 
amplitudes are typically 1.5-2 times solar. This disagreement can be partly 
explained by the lifetimes of the modes. Indeed, the oscillation modes have 
finite lifetimes, because they are continuously damped. Thus, if the star is 
observed during a time longer than the lifetimes of the modes, the signal is 
weakened due to the damping of the modes and to their re-excitation with a 
random phase. 



4 Conclusion 

The radial velocity measurements of (5 Vir, obtained over 11 nights, show a sig- 
nificant excess in the power spectrum between 0.7-2.4 mHz, centered around 
1.4 mHz, with a maximal peak amplitude of 45 cms -1 , revealing solar-like 
oscillations. Moreover, we presented the identification of thirty-one individual 
frequencies. Note that this identification is a little complicated by the presence 
of rotational splitting, as f3 Vir rotates rapidly enough to see this effect. The 
signal-to-noise of our data is unfortunately not high enough to unambiguously 
determine this splitting. The oscillation modes present an average large spac- 
ing of 72.1 yuHz, and the small spacing appears to decrease from 7.4 to 4.9 /iHz 
between 1260 and 1700 /iHz, with an average of 6.3 /iHz. As the large spacing 



is related to yM/R 3 (Kjeldsen & Bedding 1995), the mass of f3 Vir should 
be close to 1.3 M . Detailed theoretical models of (3 Vir will be reported in a 
subsequent paper. 
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